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1. Introduction 

Recent progress of computer algorithms and high performance computers enabled us to per- 
form realistic simulations of 2 + 1 flavor QCD with the physical quark masses [P. Since the equa- 
tion of state (EOS) is one of the most expensive calculations in lattice QCD, realistic simulations 
for the EOS has been performed only with staggered-type quarks [Q, ||, Q, |[|, whose theoretical 
basis such as locality and universality are, however, not well established. Therefore, to check the 
validity of these results it is important to compare with those obtained using theoretically sound 
lattice quarks, such as the Wilson-type quarks. To this end we are pushing forward a project to 
study the thermodynamics of 2+1 flavor QCD, using a nonperturbatively improved Wilson quark 
action coupled to a RG-improved Iwasaki gauge action [g]. 

To reduce the high computational cost with Wilson-type quarks, we have proposed the fixed- 
scale approach, in which we vary temperature T by changing temporal lattice size N t at a fixed 
lattice spacing a [[7L|||, ^]. Because we fix the coupling parameters for all T's, one zero-temperature 
simulation at this point can be used in common to renormalize observables at all T's. Together 
with other good features of the fixed-scale approach, we can reduce the computational cost for 
zero-temperature simulations, which is a big burden in the more conventional fixed-A^ approach. 
In the fixed-scale approach, the trace anomaly £ — 3p is calculated as usual at each temperature. To 
calculate the pressure non-perturbatively, we have developed the T -integration method' [{7|]: 

P f T jt £ ~3p 



with p(T Q ) «0. 

In this report we present the results of EOS, Polyakov loop, and chiral condensate in 2 + 1 
flavor QCD using the fixed-scale approach. Calculation of the beta-functions is updated from the 
previous reports []9[ [Tcj ]. Chiral condensate with Wilson-type quarks is known to have a severe di- 
vergence due to the explicit chiral symmetry breaking of the Wilson term [ jflj , |l2| ]. We renormalize 
the Polyakov loop and chiral condensate taking the advantage of the fixed-scale approach that the 
renormalization is common to all temperatures. 

The lattice set-up is given in the next section. Details of beta-functions and results of EOS are 
presented in Sec. 3 and 4. The Polyakov loop and the chiral condensate are discussed in Sec. 5 and 
6. The report is summarized in the last section. Part of the results including the details of the EOS 
calculation was published recently in [JS[] . 

2. Lattice setup 

A good feature of the fixed-scale approach is that we can use high-precision zero-temperature 
configurations on fine lattices, openly available on ILDG etc., for the renormalization of observ- 
ables at all temperatures. As the zero-temperature configurations, we choose those of a 2+1 flavor 



QCD spectrum study with improved Wilson quarks by the CP-PACS+JLQCD Collaboration [ |13| ] 
in the study. The action S = S g + S q is a combination of the RG-improved gauge action S g and the 
clover-improved Wilson quark action S q , 

S g = -P\Z c W^\x)+ £ • ( 2 "D 

I x,n>v x,p,v \ 
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Figure 1: Left: The global fit for j3 as a function of m p a. Square symbols show the simulation points by CP- 
PACS+JLQCD. Solid curves are the fit results around these simulation points with corresponding m p /m n 
and m^Jmtp. Open green square is the simulation point in this study. To avoid a too busy plot, only half 
of the data points are shown. Right: Beta functions on our LCP (m K /m p = 0.6337, m^/m* = 0.7377) as 
functions of j8. The scale is set by am p . Beta functions for K m i and k s are magnified by factor 100. 



s«= L E«&W> ( 2 - 2 ) 

f=u,d,s x,y 

D XJ = S Xj y - - YuP^d^y + (1 + YnWx^Sx-iiy} ~ Sx,yCswKf £ O^F^, 

where c$w is non-perturbatively determined as a function of /3 [[14]]. Among their simulation points, 
we choose j3 = 2.05, k m i = 0.1356 and K s = 0.1351, which is the smallest lattice spacing [a ~ 0.07 
fm with a scale set by rn = 0.5 fm] and the lightest u and d quark masses \m % jm p = 0.6337(38)]. 
The s quark mass corresponds to m^Jm^ = 0.7377(28). The lattice size is 28 3 x 56. 

Using the same coupling parameters as the zero-temperature simulation, we have generated 
finite temperature configurations on 32 3 x N t lattices with N t = 16, 14, • ■ •, 4 [Eh. Our range of N t 
corresponds to the range T = 174—696 MeV. The pseudo-critical point is expected to be N t « 14. 

3. Beta functions 

To evaluate the trace anomaly, we need the information of beta functions a(df}/da) and 
a{dKf/da) (f = ud and s). In this study, we define the lines of constant physics (LCP's) by m n /m p 
and m^JntQ at T = 0. The beta functions are determined nonperturbatively through the coupling 
parameter dependence of zero-temperature observables on a LCP. We use the data of am p , m n jm p 



and m^Jnity at 30 simulation points of the CP-PACS+JLQCD zero-temperature configurations | J13| ] 
to extract the beta functions. 

We have first tried to evaluate the beta functions by 'the inverse matrix method', i.e., fit the 
data of am p , m K /m p and m^Jm^ as functions of the coupling parameters j3, K U( i and K s , and then 
invert the matrix of the slopes to obtain the beta functions [^]. It turned out that a high statistics 
and many data points are required to suppress the errors in a{dKf/da) because these beta functions 
have much smaller magnitude than a(df}/da) and thus the error in the latter contaminates those of 
the formers through the matrix inversion procedure. 
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Figure 2: Trace anomaly (e — 3p)/T 4 , energy density e/T 4 , and pressure 3p/T 4 in 2+1 flavor QCD. The 
thin and thick vertical bars represent statistic and systematic errors, respectively. The curves are drawn by 
the Akima spline interpolation. 



To avoid the matrix inversion, we thus adopt 'the direct method', in which we fit the coupling 
parameters, jS, K ut j and K s as functions of the observables am p , m n jm p and m^/m^. Consulting 
the overall quality of the fits, we have refined the fit function and the fit range from our previous 
report [0]. We now fit with the third order polynomial M\. 

The result of the global fit for j8 around each simulation point is shown in Figjl] as a function 
of anip. The fit works well with x 2 /dof = 1.63, 1.08, and 1.69 for the fit of j8, K u d, and K s , 
respectively, where the degree of freedom is 10 for each coupling parameter. The resulting beta 
functions for our LCP are shown in Fig. [I] as functions of /3 . As a variable to set the scale, we 
may alternatively adopt am % , am^ or am^* instead of am p . We find that the results are roughly 
consistent with each other. Taking the results from am p as the central value, we obtain 

a^- = -0.279(24)0, = 0.00123(41)(t|), = 0.00046(26) £%) (3.1) 

at our simulation point, where the first brackets are for statistic errors, and the second brackets are 
for systematic errors estimated by the variation of the scale setting. 



4. Equation of state 

Using the results of the beta functions, we calculate the trace anomaly (e — 3p)/T 4 . See 
Ref.||] for the details of the calculation. In Fig. ||, we show the results of the trace anomaly by a 
thick red curve. 



Carrying out the T-integration of the trace anomaly according to Eq. (1.1), we obtain the 
pressure p shown in Fig. ||. The energy density e is then calculated from p and e — 3 p. Although 
our errors are still large, our EOS shown in Fig. || is roughly consistent with recent results with 
highly improved staggered quarks near the physical point [f| Eh. 
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Figure 3: Renormalized Polyakov loop and its susceptibility as functions of T, obtained with the fixed-scale 
approach. The left panel shows the results for 2+1 flavor QCD. The same quantities obtained in quenched 
QCD [ |20| ] are shown in the right panel. 



5. Renormalized Polyakov loop 

The Polyakov loop L = (3A^ S 3 ) -1 LjTrJT^ Ua{x) measures the quark free energy F(T). To 
renormalize F(T), we have to add an additive renormalization constant, which is independent of T 
and thus common for all T's in the fixed-scale approach. However, renormalization of L depends 
on T due to the relation (L) ~ e~ F l T . We adopt the renormalization scheme of Ref.[|l7|]. The 
renormalized Polyakov loop is then given as L ren = Z^ n L, where we set Z ren = 1.4801(90) from 



our potential data obtained at T = [ 1 8 ]. 

Our results for (L ren ) are plotted in the left panel of Fig. ||. As discussed in [|T0|], the In- 
dependence in these quantities are largely influenced by the renormalization factor. We note that 
our (L ren ) agree well with a result from the p4 staggered quark action in the fixed-A', approach at 
N t = 8 Ji9|], as shown by shaded symbols in the figure, although our quarks are heavier than theirs. 

In the same figure, the results for the susceptibility % Lxa are also shown. We expect the pseudo- 
critical temperature around 200 MeV. We do not see a clear peak in %l, although a peak at T ~ 1 80- 
200 MeV is not excluded due to the large errors. For comparison, we show, the results for the 
same quantities obtained in quenched QCD adopting the fixed scale approach in the right panel 
of Fig. | @. In this case, we confirm a clear peak of the susceptibility at the weakly 1st order 
deconfinement transition temperature. We thus suspect that a peak is hidden in the large errors 
and/or between the two data points. Some reasons of the lower resolution of T = {N t a)~ l in the 
full QCD study are that ./V, is restricted to be even due to the CPS simulation code we used, and 
that the lattice spacing a is coarser than that in the quenched study. We also note that the crossover 
may be weak due to the large quark mass, large N t , small aspect ratio N s /N t , etc. We reserve 
investigation of these points for future study at lower quark masses. 



6. Chiral condensate and its susceptibility 

Finally, let us study the chiral condensate. Because the chiral symmetry is explicitly broken 
with Wilson-type quarks, chiral condensate Or requires additive and multiplicative renormaliza- 
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Figure 4: Bare chiral condensates (left) and their susceptibilities (right) as a function of temperature. 



tions [11, 12], 



{W )(T) = l(TrD-i) = ^ 



(6.1) 



Here Cyy and Z ¥¥ are the renormalization constants. Because the renormalization is done such 
that the chiral symmetry is recovered in the continuum limit, we can alternatively define the chiral 



condensate in terms of the Ward-Takahashi identity [11], 

<<*(r)> 



Zpp 



(6.2) 



where P(x) and m q are the pseudoscalar density and the quark mass. Corresponding susceptibilities 
are given by 



where C,7 



Xx = C x «a#(r) 2 > - (a|(r)) 2 ) , X = wor PP, 

*2 



(6.3) 



-yy — l/Z^y and Cpp = \j (in q Zpp) L ■ Note that the factors Zyy, Zpp, and m q are just 
constants common to all T's in the fixed-scale approach. Note also that the additive renormalization 
constant Cyy cancels out in the calculation of Xw- 

In Fig.||, the bare chiral condensates and their susceptibilities are plotted. Ignoring the differ- 
ence in the overall scale, the chiral condensates and the susceptibilities are consistent with each 
other within the statistical errors. However, a clear peak could not be seen in the susceptibilities, 
as in the case of the Polyakov loop susceptibility, except for a small peak in the Xyy- 



7. Summary 

Adopting the fixed-scale approach, we have calculated the EOS in 2 + 1 flavor QCD with 
non-perturbatively improved Wilson-type quark action combined with the RG-improved Iwasaki 
gauge action. To our knowledge, this is the first result of the 2 + 1 flavor EOS with Wilson- 
type quarks. Although the light quark masses are heavier than their physical values yet, our EOS 
looks roughly consistent with recent results with highly improved staggered quarks [|3|, |5]]. We also 
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studied the Polyakov loop, chiral condensate, and their susceptibilities. With the current precision 
of the simulation, however, we could not find a clear signal of the pseudo critical point. To improve 
the situation, it is important to achieve a higher resolution in T, e.g., by performing simulations at 
odd N t and/or by combining data at different scales. We are planning to test these ideas in a study 
at the physical point, using the on-the-physical-point configurations generated by the PACS-CS 
Collaboration [p. 
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